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Abstract 



We discuss the concept of connected, reparameterization invariant matter correlators 
in quantum gravity. We analyze the effect of discretization in two solvable cases : 
branched polymers and two-dimensional simplicial gravity. In both cases the naively 
defined connected correlators for a fixed volume display an anomalous behavior, 
which could be interpreted as a long-range order. We suggest that this is in fact 
only a highly non-trivial finite-size effect and propose an improved definition of the 
connected correlator, which reduces the effect. Using this definition we illustrate 
the appearance of a long-range spin order in the Ising model on a two-dimensional 
random lattice in an external magnetic field H, when H —>■ and (3 = (5c- 
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1 Introduction 



In a theory where gravity is quantized it is non-trivial to define the concept of a con- 
nected correlation function. The problem is only apparent once one has a genuine 
non-perturbative definition of quantum gravity where one can go beyond the expan- 
sion around flat space. To exemplify the problem let us define a reparameterization 
invariant correlation function in quantum gravity: 

G^\R) = J V[g) J V[X) e~ SA J J dV 9 {t)dV g {0 MOMO 6{D g (t, £') - R). 

(1) 

Here dV g is the invariant volume element and -D 9 (£, £') the geodesic distance between 
£ and £' for a given metric g. T>[X] symbolizes the integration over additional degrees 
of freedom of the theory, i.e. the matter fields and (pi are local observables built from 
these fields and/or the gravitational field. The functional integral is over equivalence 
classes of metrics, i.e. the group of diffeomorphisms is divided out. Finally A is the 
cosmological constant, i.e. we have the following decomposition of the action: 

S A = AV g + S v [g,X], V g = J dV g {Z), (2) 

where Sy is independent of the cosmological constant. As we can see from (|IJ) the 
correlator involves three rather than two non-trivial operators, contrary to the flat 
space case, where the definition of a distance does not involve an extra operator. 
The new element is the geometric "separator", which for any non-trivial geometry 
becomes a complicated non-local object. A partition function for quantum gravity 
can be defined by 

Za = J V[g\ J V[X\ e~ s \ (3) 

and Gf 02 (i?) is an unnormalized correlation function. The simplest object of this 
kind is a "geometric" two-point function, including only the "separator" : 

G\\R) = J V[g) J V[X) e"** J J dV g {i)dV g {i') 5{D g {^') - R), (4) 

i.e. the same object as in ([I]), just with 0(£) — > 1. Gj^(R) can be viewed as the par- 
tition function for the ensemble of universes where two marked points are separated 
by a geodesic distance R. We have 

j o d RG\ 1 (R) = ^Z A = Z A (5) 

This two-point function plays a special role in describing the geometric properties 
of the system 0. In general Gj^lR) is expected to fall off exponentially, reflecting 
the fact that there is an exponentially small probability to create a universe where 
two points are separated by a geodesic distance R much larger than some power of 
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the cosmological constant. This power sets the geometric scale of the systemfj, or a 
relation between the average radius (R) and the average volume (V) of the universe: 
(V) oc (R) dh , where d h is the Hausdorff dimension. The concept of a geometric scale 
(R) is very important when we discuss the problem of matter correlations, which in 
principle may involve some other "physical" scale. The geometric scale is controlled 
by the cosmological constant, while the physical scale will in general depend on some 
other coupling constants, describing the matter sector of the theory. Away from the 
critical point in this sector we expect the physical scale to be small compared with 
the geometric scale. Only then we can expect to distinguish those two scales In the 
following we shall always assume that it is the case. With a small abuse of notation, 
we call this limit the thermodynamic limit. Close to the critical point the two scales 
become comparable and there we may expect to find a non-trivial coupling between 
the matter and gravitational sectors of the theory. If we view Gj^ (R) as a partition 
function, the distance R plays the role of an additional coupling constant of the 
theory. All these peculiarities make the definition of connected correlators difficult 
and sometimes ambiguous. 

A priori we have various possibilities of defining a normalized correlation func- 
tion. Let us mention two different definitions: 

(0 X 2 ) A (R) = ^V^Z (6) 



Z 



A 



and 



<fcfc>A(«) = -^p (7) 

The definitions differ by the way geometry is counted. In (0) it is counted in the 
same way in numerator and denominator. In fact, if <pi and 2 are independent 
observables which do not couple to gravity at all, definition (H) gives 

(0 1 2 ) A ( J R) = (0 1 )(0 2 )^1, (8) 

On the other hand, the definition (^) leads to 

(0102> A (i?) = <01><02), (9) 

independent of R. 

The results (§) or (||]) are trivial from the point of view of the correlation func- 
tions. Both should be viewed as examples of a situation, where the observables <pi 
at different points are uncorrelated. The correlator can in this case be expressed as 



4 The definitions (@)-(|) apply to any field theory of quantum gravity. However, when we say that 
the geometric scale is set by the cosmological constant, we have in mind two-dimensional quantum 
gravity. In four-dimensional quantum gravity the geometric structure might be more complicated 
since the gravitational coupling constant is expected to set the scale of quantum fluctuations and 
the cosmological coupling constant the scale of the universe. 
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a product of three single-operator averages, which is best seen in (H). One would 
expect a similar structure of the correlation function in the thermodynamic limit 
when the observables 0j do couple to gravity, but when the distance between the 
two points R is much bigger than the correlation length. The expectation values 
(4>i) should then be replaced by (<fti) A implying the non-trivial dependence on the 
cosmological constant. The R dependence should factorize, like in @ or (||). 

At smaller distances there may be deviations from a simple factorization []. The 
obvious questions are: Can these be interpreted as a signal of a correlation between 
the fields? Which type of behavior can still be attributed to the uncorrelated opera- 
tors? Such questions clearly can not be answered without some hint from a solvable 
model. Unfortunately the number of solvable theories, where one can actually com- 
pute the two-point functions, is very limited. In the next section we shall present 
some exact results concerning operators which we believe are uncorrelated. 

Of course a behavior like (|8|) and (0) can not be observed if we have correlators 
with infinite correlation length (or where the geometric and physical scales become 
comparable), as will be the case when we consider correlators between the conformal 
fields coupled to two-dimensional quantum gravity. 

To obtain a consistent definition of the matter correlation length we must define 
the concept of a connected correlation function. Various definitions have been given 
in the literature || ^, |5| and below we summarize the discussion. First we face the 
problem of a sensible definition of (</>) A . Several choices seem possible, but in order 
to match the definitions of correlation functions, as given by (||) or (0) one can use 
either 

= tJ v[s] I v[x] e " SA "./"wwfl - §§wm (10) 

or ^ 

(0> A (fl) = ^l) A (i2) = ^^ (11) 

The last definition depends on the geodesic distance R, and it is from this point of 
view slightly unusual as a definition of an expectation value of a field. But not much 
more than in usual field theory where the lack of translational invariance caused by 
an external field can introduce a space-time dependence in ((f)). As discussed above, 
such dependence is expected to disappear for large R implying the existence of a 
well defined limit {4>) A (R — > oo). In fact, this limit should agree with the value 
defined by (|ToD provided the correlation length of the matter fields is (much) less 
than the average radius of the universe. The definition of (0) A is related to the 
correlation between the and unit operators, as is clear from (|H]) and (0). Thus 
the statement that the two definitions of (0) agree in the thermodynamic limit, for 
large R is equivalent to the statement that we have a factorization 

Gl\R) **(<!>) A G%(R) (12) 



3 An attempt to derive a systematic operator product expansion can be found in reference 
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for such values of R. In general we would expect the unit operator to be uncorrelated 
with any other operator, so naively we would expect ( JP2| ) to be satisfied almost for 
any value of R. This can be checked numerically. 

Note that the correlation function G'^(R) almost inevitably enters in any sensible 
definition of the connected part of a correlator since one has to consider an object 
like 

(0i0 2 )7 (R) = ((0i " (0i> 1)(02 - (02) 1)) A (R), (13) 

with some definition of (0,). If we use ([IT]) as a definition of (0j) and (|7|) as a 
definition of the correlator, then eq. ( |i~3D fulfills the standard decomposition and 
can be written as 

(0 1 2 ) A ( J R)-(0) A (i?)(0 2 ) L (i?). (14) 
If we use (H) and ( pTj|) instead, eq. (|Tj|) can be written as 

(r) = (is) 

i- (Gt^(R) - (<h) A Gp(R) - (fa) A Gfi(R) + (0 X ) A (0 2 ) A Gi 1 (i?). 

In this case one does not have the standard local decomposition as in eq. (0), but 
integrating over R one obtains: 



) A -(0l) A (0 2 ) A , (16) 



where the first term in (p7j) is the usual integrated correlator as known from two- 
dimensional quantum gravity: 

= * J V[g] J V[X] e~ s * J J dV a (0dV g (O 0(00(0- (17) 



We can now define the correlation length by the exponential decay of (00) « (R), 
defined either by (|i~4] ) or (|15|) and the thermodynamic limit is where this correlation 
length is much smaller than the average radius of the the universe. 

The fixed volume partition functions Zy and Gp^ 2 (R) are related to Z\ and 
G' A ^ 2 {R) by Laplace transformations 



dVe~ AV Z V: Gl W2 (R)= dVGl l92 (R). (18) 

Jo 

and it is natural to use the definition of correlators, expectation values of fields and 
connected correlators corresponding to eqs. (H), (§§), (|T0|) and ([15]), just with the the 
partition functions for fixed A replaced by the ones for fixed volume V, given by 
([18]). Note in particular that we have 

(0)y=^/W(/^ 5 (O-^) jv[X]e- s ^^jdV g {t,)<j>{X-i), (19) 



5 



as one would have expected. The connected correlator for finite volume could thus 
be defined as 



(^2)7 (R) = (20) 
— L. (Gp^(R) - (<h) v G»*(R) - (&> v Gj*(J2) + (<h) v (<h) v GvW 
with the normalized fixed volume correlators defined as 

G^(R) = y^G^R), (21) 

satisfying 

roo _ 

dRG l v \R) = V. (22) 



As we show in the next Section, the discrete regularization of a theory may, 
and in fact does lead to some complications, where the finite volume effects tend to 
mimic the physical correlations even in situations where there are no correlations in 
the grand canonical formulation. This was first realized in [||. 



2 Analytical results 

There are unfortunately very few systems, where the concepts presented above can 
be compared with the analytic prediction. We present here the few models where 
the two-point function can be explicitly calculated. These systems are branched 
polymers |6|, |7| and two-dimensional simplicial gravity |{|. In both cases one 
defines a discretized (integer) geodesic distance r between the points of the manifold. 
The physical relation between the the continuum and discrete geodesic distances is 

R = ra, (23) 

with the lattice spacing a — > in the continuum limit, but the distance R kept fixed. 
In both cases the relation between the the continuum volume V and the discrete 
volume N is anomalous Q: 

V = Na d \ (24) 

where dh is the Hausdorff dimension of the system, equal 2 for a generic branched 
polymer and 4 for two-dimensional simplicial gravity. 

In both cases there is no extra matter content in the theory and the only observ- 
ables we can discuss are related to the local geometric properties of the manifold. 
In the case of simplicial quantum gravity these can be some functions of the coor- 
dination number of a vertex, in the case of branched polymers - functions of the 
branching ratio at a vertex. Both cases correspond to the observables, which we 
expect to be essentially uncorrelated for r > 0. 

6 One usually chooses the scaling R = a 2 / dh r and V = a 2 N rather than ( p3| ) and (p4|). However, 
( |25| ) is more convenient from a notational point of view in the arguments to follow, so we prefer 
to work with a "rescaled" cut-off a defined by ( p3| ) and (p4|). 
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2.1 Branched polymers 

Let us start by repeating the discussion of the simpler case of the branched polymers 
H [/J. The partition function in this model is given as a weighted sum over the 
ensemble of trees. Trees are weighted by one- vertex branching weights p(qi). The 
partition function for the ensemble of trees is given by : 

Z = £ exp(- M JV) E T^Y IT Pin) (25) 

where is the order of a vertex z, iV is the number of vertices and C(T) is an 
appropriate symmetry factor of the graph; \x plays the role of the bare cosmological 
constant. 

The correlation functions are constructed by means of the partition function 
Z of planted, rooted, planar trees. This partition function can be found from the 
following recursive equation || : 

Z = e-»F(Z) (26) 

where 

F(Z) = J2p(q)Z«-\ (27) 

q=l 

In the generic case equation fl26|) has a critical point \ic for which 

e-T(Z ) = 1, Z = Z(nc). 

For /i approaching the critical value \xc from above, Z has the following singularity : 

Z([i) = Z -Z' 1 ^jI=jr c + O{v-pi c ) (28) 

The natural definition of a distance r between two vertices of a graph is the 
number of links joining these vertices. The discrete analogue of the geometric two- 
point function can be calculated in terms of Z(p) 0, 0] as 

<w = (i + -§^V z% - ( 29 ) 

The factor Z 2 in ( p9[) is a contribution from the two ends of the chain. Close to the 
critical point 

A = -log(l + — ) oc ^/jT=JT c + 0(n- n c ). (30) 

and ( p9[) becomes 

Gi 1 (r) = exp(-(r-l)A)Z(A) 2 . (31) 

In a branched polymer system the only local observables we can construct are 
functions of vertex orders q\, q2 at the end points. Replacing Z — > Z e Xi F(Ze Xi ) / F(Z) 
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at the end points we obtain a generating function of all such observables: differen- 
tiating with respect to A, at A, = we generate powers of % = 1,2. Notice that 
for Z satisfying fl26|) the generating function 

where the averages are taken with respect to the partition function Z{\x). We con- 
clude that for every choice of observables we have a simple factorization 

G^(r)=g(X 1 )g(X 2 )Gf(r), (33) 

which proves a lack of correlation between any pair of vertex order operators at 
points 1 and 2. 

Near the critical point the generating function (7(A) can be expanded as 

(7(A) = 0o(A) + A&(A) + 0(A 2 ) (34) 
Using the simple form fl3"T| ) of G^(r) we get 



Gj*.M( r ) = ^ (A 1 ) + ^ 1 (A 1 )a r + ---)(^o(A 2 ) + ^i(A 2 )9 r + ---)G; i ( 



^o(Ai)^ (A 2 )Gi 1 (r + 5(A0 + <J(A 2 )), (35) 



where we have introduced a s/w/t 



MA) = (36) 



and where the dots stand for higher order terms typically proportional to d 2 G]} (r) / dr 2 . 



2.2 2d gravity 

Two-dimensional simplicial gravity can be obtained as the planar limit of the <p 3 
matrix theory The coupling constant g = of this theory can be parametrized 

as 

8g 2 = s(l-s 2 ). (37) 

In this parameterization g = corresponds to s — 1 and the critical value g 2 c = 
l/12v^3. In the planar limit each 3 vertex is dual to a triangle and a 3 graph can 
be viewed as a two-dimensional surface built from triangles, expansion in powers of 
g becoming the expansion in the area of the surface. Following || we construct the 
transfer matrix [10| 



GpiX, Y;r) = J2 X L Y L 'G,(L, L'\ r), (38) 



L,L' 
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where G^L, L'; r) is the sum of all possible connected planar 3 graphs with bound- 
aries being the loops (planarly ordered sequences of external links) with lengths L 
and L', separated by a "distance" r. The distance between the two vertices of a 
graph is the length of the shortest path following the links of a graph. The transfer 
matrix G^(X, Y; r) can be calculated using the "peeling" method JTI] giving 

G > {m = WT)^ (39) 

where 




and X(X,r) satisfies 



/„(X, 9 ) = (1-— (40) 



4logS4 (41) 



A & T(X) 




\A 2 ~ A 9X 



A 

T{X) 

For g — > g c = e~ Mc we have 

[i — fx c oc A 4 + • • • 

which reflects the fact that the Hausdorff dimension dh = 4. Introducing T = T(X) 
and T(X) = e~ Ar °W we have 

f _ e -A(r+r (X))_ ^2) 

Using this parameterization we easily express f^{X,g) as 

f»(X,g) = ^coth(jA(r + r ))(l-coth 2 (jA(r + r ))) (43) 
g 2 2 2 

= F(A(r + r )). 

The transfer matrix described above can be used to calculate the two-point correla- 
tors 0, §. Below we shall discuss only the simpler case, when one of the operators is 
the unit operator. The solution for the general case with two non-trivial operators 
can be deduced from symmetry. The first step will be to close the incoming loop. 
We are still left with the open final loop. The resulting function will be denoted by 
G^X; r) and can be expressed as 

dY C (X l -rW 9 



G^X-r) = t 7 ^ 7 G M (X,-;r)Y— *(y) (44) 



jjrfM^' (45) 
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where 

IX . 1 sX, I AgX 

»W = j( 7 + D + i(i- 7 f-^ («) 

is a generating function of the connected Green's functions of the 3 matrix model 



[16|. For small A the effect of closing the incoming loop can be represented as 



f, , Y \ F(A(r + r ))-d r F(A(r + r )) + --- 

G ^ X]T) ~ FWo) (47) 

F(A(r-l+r )) 
F(Ar ) 

The dots correspond to terms proportional to the second derivative of F. The 
X dependence in this formula appears through tq = ro(X). This is exactly the 
information we need. The particular choice of a set of operators we wish to study is 
not very important, provided the incoming lines represented by the X dependence 
are attached to some local geometric object. The simplest choice is just a line, 
joining the end points, which leads to the generating function 

G^(r) = G (1 (X = je Al ;r) ) (48) 

where we introduced, as before, the parameter Ai to count the length of the line and 
its higher moments. Notice that although the r dependence of this quantity is more 
complicated than in the branched polymer case, the fundamental properties remain 
the same. Denoting by S(X) = r (A) — r (0) we have 



Gl X ^ =0} (r)=g(\ 1 )Gf(r + 5(\ 1 )). (49) 

F(Ar (0)) 
F(Ar (A))' 

Again the correlation function involving the nontrivial operator (A 7^ 0) is related 
to that of the unit operator (A = 0) by a A-dependent multiplicative factor and a 
A-dependent shift of r. 



where 

0(A) 



2.3 The Ansatz for a connected correlator 

There are two important lessons one can learn from the examples presented above. 
The first one is that even in cases where there is no correlation the naive factorization 
of the two-point correlator: 

GlHr) = (<!>}, Gf(r) (50) 

is satisfied only asymptotically, for large enough r. This does not necessarily mean 
that there is no simple factorization, as we see from the example of pure 2d gravity, 
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(see (42)). The object for which we observe a simple behavior is not the two- 
point function itself, but rather a geometric object, natural to the evolution of 
the characteristic equation of the transfer matrix. The two-point function is itself 
expressible in terms of this object. Near the critical point the geometric scale of the 
system is provided by A and the relation between the cosmological constant and 
this parameter is anomalous 

H - n c oc A dh . 

In general only for r large enough the correlation function Gj^(r) has a simple 
exponential form. 

The second lesson is that we can nevertheless extend the idea of factorization to 
be satisfied for the whole range of r if instead of (50) we use the modified Ansatz: 

G^(r) = (A,~B^d r + ---)G l l !(r) (51) 

which asymptotically (for r large enough for Gj}(r) to be purely exponential) cor- 
responds to 

= ^ + A + • • ■ 
To order 0(A) the Ansatz (|5lD can be viewed as an additional shift: 

G i ;{r) = A 4> G 1 ^r + 5 !t> ) (52) 

with 

**- V 

For two uncorrelated operators the corresponding Ansatz would be 

G^{r) = Afa A^ Gf(r + 5 1 + 5 2 ) + 0(A 2 ). (53) 

The additional shifts are specific for the operator and are additive. We also made 
an assumption that A^. ^ 0. The uncorrelated contribution described above has to 
be subtracted if one is interested in the connected correlation function. The effect 
described above is a finite-size correction. In the continuum limit r — > R/a and 
d r — > a Or. It is however very important if we analyze the system using the fixed 
volume canonical numerical simulations. 

The unnormalized correlation function G^ 1 ^ 2 (r) can be represented as a (discrete) 
Laplace transform of the fixed- contributions as 

GfHr) = EexpC-zJV) C N G**»(r), (54) 

N 

where Cn controls the normalization of the fixed- volume correlation function and is 
independent of 0. The commonly used normalization corresponds to 

£G#(r) = JV. (55) 
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Using standard arguments Gp (r) can be obtained by inverting the Laplace trans- 
form (EQ). From (EBi) it follows that 



GT 2 (r) = [Afa -B^d r + - ■■)(A^-B^d r + - • -)G%(r) w G 1 1 ^(r + 6 1 + 5 2 ). 

(56) 

In the large- N limit we expect 

G#(r) « tf-Mgix), (57) 

where the scaling variable x = (r+8o)/N 1 ' dh . It is clear from ( [56] ) that derivatives d r 
produce finite-size corrections 0(1/N l l dh ) or 0(1/ L), where L is the linear (discrete) 
size of the system. 

In a numerical experiment at a fixed volume N one measures the correlation 
functions G^(r), G]f l (r) and G^ 2 (r), where G^(r) is normalized as in (|55|). As 
we argued above we expect to 0(1/ L 2 ) 

G)} i (r)=A <j>i G%(r + 5 i ). (58) 

This behavior can easily be checked in numerical experiments and we find it satisfied 
with amazing accuracy, not only in cases discussed above, but also for other two- 
dimensional systems, as will be discussed in the next section as well as in correlation 



experiments in four-dimensional simplicial gravity [17|. Following (|58"D one would 
expect Ai to be an iV-independent constant A^. = (<pi) N=00 - Experimentally we 
found that an excellent fit to ([58]) corresponds to the choice A^. = (4>i) N , which for 
an extensive quantity differs by 0(1/N). 

In general the corrrelation function G^ 2 (r) may contain the non-trivial con- 
nected part. However if we define the connected part of the correlator in the analogy 
with (gOD as 

G%^ 2 ' conn (r) = (59) 
GftHr) ~ (<Pi) N G]f 2 (r) - (th) N G]t(r) + 
(h} N (h) N G^(r) } 

even in case when there is no correlation in the grand canonical ensemble, as dis- 
cussed above, we find a non-zero contribution 



% d 2 r G%(r) = B 1 B 2 N 1 - 3 ^f(x), 

where f(x) = g'(x), x = (r + 5o)/N l l dh . The measured correlation mimics a real, 
physical scaling correlation function, corresponding to an object with the anomalous 
dimension 2/dh- Note that this behavior will be observed even if the two-point 



Following the discussion above it represents in fact only the contribution from the 



function Gj/(r) is purely exponential, as is the case for the branched polymer system. 
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disconnected part of the correlation function and can be viewed as an artifact of the 
Laplace transform to a finite volume N. It is exactly this part, which we would like 
to eliminate, when we try to measure the connected part of the correlation function. 

Our discussion shows that it may be impossible to eliminate completely the 
disconnected contribution in a finite volume correlation experiment. However, we 
may reduce it. The following simple redefinition of the connected correlator: 

G% lH ' conn (r) = (61) 

(h) N (4>2) N G^(r + 6 1 + 5 2 ), 

gives an extra factor N~ 2 ^ dh as compared with (|5^) and in all practical cases was 
sufficient to reduce the disconnected signal below the error level. This point will be 
discussed further in the next section. 



3 A numerical recipe 

In this section we present a simple algorithm, which permits us to measure non- 
trivial correlators in numerical experiments. Assume that from the experiment we 
know four correlation functions : G^ 2 (r), (r), (jjy (r) and G]}(r). To make 
use of the formula (fOP we need four parameters: (<fii) N , (^2)^, &i and <5 2 - 
As a first step we use 



G^(r) = (4 >l ) N G^(r + 6 l ). (62) 

Correlation functions measured in the experiment are given only for integer values of 
the distance r. To interpolate between these values we use a four-point interpolation 
formula, which is used to obtain a function for k > r > k + 1 from it's values at 
k — 1, k, k + 1 and k + 2. Two parameters (4>i) N and 5i are fitted by the least 
square method. We decided to use this method rather than the y 2 since including 
the numerical errors has an effect of giving a large importance to the ends of the 
distribution, where both the correlation functions and errors are small. In practice 
it turns out that the fitting is fairly insensitive to the extrapolation method (we 
used the spline method as an alternative). The value of (4>i) N obtained in a fit was 
practically the same as the average obtained independently from the experiment. 
After some checks we decided to use the experimental value from the start and to 
use (B2p only for a one-parameter fit to find the optimal value for One should 



stress that in all cases the fit is very good provided (<f>i) N is not to small. From 
our discussion it follows that the connected part of the correlation function remains 
unchanged if we add a constant to the operator. We used this property, when 
necessary to avoid numerical difficulties. 

We start our presentation with the numerical results for the branched polymer 
model for a system with 1024 vertices and weights p(q) = 1/q. The results for 
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Figure 1: Branched polymer. 



the connected correlators between the branching ratios at two points separated by 
a distance r defined using ( |59| ) and ([H]) are presented on the figure [1| This ex- 
ample illustrates the power of the proposed algorithm. We can see that a strong 
disconnected signal is reduced almost completely and that we are left with only local 
correlation at r = 0. 

Similar effect is presented on the figure |2| for pure 2d gravity using the data 
from the combinatorial spherical surfaces with regular triangulations and with 8000 
points. The definition of a distance used in this measurement is different than in the 
discussion presented in the section 2. Instead of the distance between the 3 vertices 
(centers of the triangles of the surface), defined in terms of dual links, we use the 
distance between the vertices of the surface, defined as a shortest path following the 
links of a surface. We measure the connected correlator of the orders of two vertices 
(numbers of triangles containing this vertex). This example shows that also in this 
case we observe the anomalous long-range behavior if the the connected correlator 
is defined using (^). Using an improved definition (|UT]) we observe a nontrivial 
correlation only for very small r. 

We used the same method to study the spin-spin correlations for the ferromag- 
netic Ising model in a non-zero magnetic field H on a random lattice. This is an 
example of a theory with a non-trivial matter content. The action for the spin sector 
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Figure 2: 2d gravity. 



of the theory is chosen as 



ft = -E(E^,.i + E ff ( 63 ) 

The spins are placed in the centers of triangles and we use degenerate triangula- 
tions which allow the two vertices to be connected by more than one link and the 
links connecting a vertex to itself. This model was solved in [[0J. For H > the 



model is always in the ordered phase with the geometric properties of pure grav- 
ity (jstr = —1/2 and d h = 4). For H = the model has two phases depending 
on the value of (3. For (3 = j3 c = — log((2v7 — l)/27) it undergoes a third order 
phase transition between an ordered and a disordered phase. At the transition the 
geometric properties change (jstr — — 1/3). In our experiment we measure the cor- 
relation between 0j = (1 + Sj)/2, (which is 1 for up spin and zero for the down spin). 
This choice avoids problems close to H = 0, where (sj) — > for every f3 on a finite 
lattice. In our case — > 1/2. The connected correlator is independent on the 
additive constant and is (up to a trivial factor 1/4) simply the spin-spin correlation 
function. The distance used is the triangle distance, i.e. the shortest path on a dual 
lattice. We can predict that if the algorithm works, the observed shift 5(H) should 
vanish both for large H and for H — > 0. For large H the system becomes completely 
ordered and the correlation function Gff(r) becomes equal to G]j(r). For H — >• 
average spin approaches zero and G]f(r) ~ \@n(t)- Below we show results from 
numerical experiments for system sizes 2k, 4k, 8k and 16k triangles performed at 
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Figure 3: Spin-spin correlation functions at H = 0.0125 defined by fl59l). 

(3 — (3 C for various values of the magnetic field H. Similar extensive experiments at 
H = were performed by ||18|| , although the definition of distance was different. 

On figure ^|we present the connected correlation functions Gff(r) versus r for a 
"large" magnetic field H = 0.0125 defined in the standard way fl5"9"|). On this and 
consecutive plots we use colors to distinguish the system size: black for 2k, red for 
4k, green for 8k and blue for 16k. We see a large disconnected contribution, which 
in fact becomes stronger with increasing the system size. 

On the next plot we show the same functions scaled by a factor N z / dh ~ l with 
dh = 4 plotted versus the scaling variable x = (r + 5o)/N l / dh . 

For the sake of presentation we decided to use the theoretical value = 4 
rather than to fit this value. The shift 5q was obtained by matching the scaling of 
the G]}(r) correlators, again assuming dh = 4. As we can see the long-range part 
of the correlator scales exactly as predicted by fl5D|). After reduction by ( |5"T| ) the 
correlators contain only the non-scaling part. We show it in figure |5|. The shift 
6(H) = .35 ± .02 is for this value of H practically independent of the system size. 

Decreasing the magnetic field H brings us closer to the phase transition. When 
we apply the algorithm presented above we discover that for each volume the shifts 
first grow, and then decrease back to zero. The value of the shift as well as the 
position of the turn over point start to depend on a system volume. From the scaling 
arguments we expect that near the critical point we should observe a universal 
behavior of spin-spin correlations, provided the magnetic field is taken as volume 
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Figure 6: Shift 6(h) vs. h. 
dependent in the following way]] 
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On figure | we plot shifts <5at(/i) for the four system sizes. We see that the position 
of the turn-over point corresponds roughly to h = 1. Dashed lines join points with 
the same value of H. 

Below h = 1 we observe also a change of 7^ from the pure 2d gravity value —1/2 
to the c = 1/2 value 7 str = —1/3. We measured this parameter using the standard 



method of measuring the minbu distribution [[L9| poj] . 

Finally figure [5] represents the correlation functions Q^' conn obtained using the 
improved formula (fHf), scaled by a factor jV 1 / c( ' i+1 / 3_1 and plotted versus the scaling 
variable x = (r + 5o)/N l l dh . Values of dh = 4 and So = 4.5 are the same as for 
other plots. Plots correspond from top to bottom to different values of the scaled 

7 In flat space we expect from standard scaling arguments the singular part of the free energy 
behaves as 

F sm (0, H) = F sm (l,H/9^ d+ ^/ 2 ), 9 = \(3 b ~ 0\, 

where j m is the critical exponent for the magnetic field. This formula is converted into a finite 
size scaling relation by using that the pseudo-critical point for volume V is obtained when the 
correlation length £ = Q~ v is equal to the linear dimension of the system, i.e. to L = V 1 l d . Thus 

F sm (V,H) =F sm (l,HV^ 1+ ^). 

Thus we expect a universal dependence of H Vi( 1+: %r\ The exponent jm/vd —2/3 after coupling 
to quantum gravity. 
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Figure 7: Parameter 7 str vs. h. 



magnetic field h = 0.2, 0.4, . . . 2.0. The distributions were obtained in practice by 
interpolating the measured correlation functions at weak fields H to match the scaled 
values h. Indeed we observe scaling as predicted. 



4 Discussion 

A reasonable definition of a connected correlator is given by eq. for a fixed 
cosmo logical constant, and by eq. (|20| ) for a fixed space-time volume. However, 
even when no correlations exist for fixed cosmological constant, the naive discretized 
analogue of (^), as defined by eq. ([59]), has a non-trivial scaling. This is made 
explicit in eq. (^). Basically a dominant disconnected part is still present in the 
definition (|59|). As we showed, the refined definition (^1|) manages to cancel this 
dominant disconnected part and leaves the genuine connected part of the correlator 
as the dominant part of the two-point correlation. 
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Figure 8: Extrapolated spin-spin correlation functions for (from top to bottom) 
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